Transfer matrix method for optics in graphene layers 
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A transfer matrix method is developed for optical calculations of non-interacting graphene lay- 
ers. Within the framework of this method, optical properties such as reflection, transmission and 
absorption for single-, double- and multi-layer graphene are studied. We also apply the method to 
structures consisting of periodically arranged graphene layers, revealing well-defined photonic band 
structures and even photonic bandgaps. Finally, we discuss graphene plasmons and introduce a 
simple way to tune the plasmon dispersion. 
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I. INTRODUCTION 

Graphene is a fiat monolayer of graphite with car- 
bon atoms closely packed in a two-dimensional honey- 
comb lattice. A hallmark of graphene is the existence of 
Dirac cones in the electronic band structure, resulting in 
extraordinary structural and electronic properties with 
great potential for nanoelectronics PHI]- 

In addition to outstanding electrical, mechanical and 
chemical properties, graphene has interesting optical re- 
sponse. One of the striking optical properties of graphene 
is that its reflectance, transmittance and absorbance are 
determined by the fine structure constant 0, Q . Despite 
being only one-atomic-layer thick with negligible reflec- 
tion, a single free-standing graphene shows significant ab- 
sorbance, universally about 2.2% in a spectral range from 
near-infrared to visible [H H, 0] • In the infrared regime, 
graphene absorption can be altered by applying gate volt- 
ages [8|,|9(. For few-layer graphene, the optical absorption 
is proportional to the number of layers [lOj , leading to a 
visual image contrast which can be used practically to 
identify the number of graphene layers on a substrate. 
The highly transparent and outstanding electrical prop- 
erties of graphene make it attractive as transparent elec- 
trodes [n], [l2j ■ The broadband absorption implies that 
graphene has the potential as an active medium used in 
broadband photodetectors fl3l Il4j , ultra- fast lasers [l5[ 
and optical modulations [16| |. 

In doped or gated graphene, collective excitations — 
plasmons exist with interesting optical features such as 
dee p su bwavelength and high confinement of optical fields 
Il7l423| | , similar to surface plasmons in metal surfaces (23 - 
|26[. As a result, graphene may serve as a one-atom- 
thick platform for infrared and terahertz metamaterials 
[27l [28j . A number of photonic devices such as waveg- 
uides, splitters and combiners and superlenses could be 
envisioned |28l . l29( . Numerical simulations suggest that 
periodically patterned arrays of doped graphene nan- 
odisks may completely absorb infrared light at certain 
resonant wavelengths [30( , which was soon confirmed ex- 
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perimentally 31( . These interesting optical properties of 
graphene may also offer potential applications in photon- 
ics |32j |. 

In this paper, we develop a transfer matrix method for 
studying optical properties in non-interacting graphene 
layers. This paper is organized as follows. In section 
ILT1 we introduce the transfer matrix method for study- 
ing the propagation of light through graphene layers, to- 
gether with the optical conductivity of graphene used in 
our calculations. In Sees. HTT]to|Vl we apply the transfer 
matrix method to the study of optics in graphene layers. 
Specifically, in section [TTTI we discuss reflection, transmis- 
sion and absorption in single-, double- and multi-layer 
graphene. In section HVl we discuss photonic band struc- 
tures in periodical graphene layers. In section [V] we 
discuss plasmons in graphene. Finally, we present our 
summary in section fVTl 



II. GENERAL FORMULATION 

Transfer matrix method is a powerful tool in the analy- 
sis of light propagations through layered dielectric media 
[33l Hi| . The central idea lies that electric or magnetic 
fields in one position can be related to those in other po- 
sitions through a transfer matrix. Within the framework 
of the transfer matrix method, there are two kinds of ma- 
trices: one is the transmission matrix that connects the 
fields across an interface and the other is the propaga- 
tion matrix that connects the fields propagating over a 
distance within a homogeneous medium. 



A. Transmission matrix 

We first consider the propagation of light across an in- 
terface formed by a graphene layer which separates two 
dielectrics of dielectric constants E\ and £2, as shown 
schematically in figure [IJa) . The graphene layer has an 
optical conductivity a lying at z = 0. Light is assumed 
to be polarized in the y direction and propagate in the 
z direction. For structures considered, s and p polariza- 
tions can be decoupled. As a result, we can deal with s 
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FIG. 1. (a) A single graphene layer surrounded by two 
dielectrics of dielectric constants ei and si. The graphene 
layer is characterized by a conductivity a. Red and purple ar- 
rows indicate incoming and outgoing light, respectively, (b) A 
stack of N graphene layers of conductivity en (i = 1, 2, N) 
are separated by different dielectrics of dielectric constants 
Si (i — 1, 2, N + 1). The spacing between two adjacent 
graphene layers is denoted by (i = 1, 2, N — 1). 



and p polarizations separately. 

For p polarization, the magnetic field is polarized along 
the y direction and can be written as the form 

H ly = (a ie tkl * z + b ie - lkl * z y kl * x , z < 0, (1) 
H 2y = (a 2 e lk ^ z + b 2 e~ lk ^ z )e lk ^ x , z > 0. (2) 

Here, a; and bi (i — 1,2) are the field coefficients, ki x 
(ki z ) is the x (z) component of the wave- vector ki = 
y/FiUj/c, where uj is the angular frequency and c is the 
speed of light in vacuum. The first (second) term in 
the parenthesis on the right side represents propagating 
waves along the z (—z) direction. From the Snell's law, 
we will immediately have k\ x = k 2x . 

The electric and magnetic fields at the interface can be 
related by the following boundary conditions [35[ 



nx (E 2 -E 1 )| 2=o = 0, 
nx (H 2 -H!)U =0 = J, 



(3) 
(4) 



where n is the unit surface normal and J is the sur- 
face current density of the graphene layer. Applying the 
above boundary conditions at z = 0, we will have 



ki z 



(ai 



ko 



(a 2 -6 2 )=0, (5) 

£l £2 
(ai + bx) - (a 2 + b 2 ) = J x . (6) 

Note that J can be obtained from Ohm's law, namely 

ak 2z 



Jx — vEx\z=o — 



s e 2 uj 



-(a 2 - b 2 ), 



(7) 



where £o is the vacuum permittivity. Combing equations 
©-([7]), the coefficients a\ and bi can be related to a 2 
and b 2 by a 2x2 transmission matrix D\^ 2 



a i 



a 2 



where 



Dl^2 = ~ 



Vp 
Vp 



1 



-tp 

<p 1 



Vp 
Vp 



(8) 



(9) 



with the parameters r\ v and £ p given by 
£\k 2z „ <rk 2z 



Vp 



e 2 k x 



£ £ 2 W 



(10) 



For s polarization, the electric field is polarized along 
the y direction. Similarly, by applying the boundary con- 
ditions and Ohm's law, the transmission matrix for s po- 
larization that relates the electric fields at the two sides 
of the interface can be obtained, as 



>2 — 7T 



i - Vs - 6 i + Vs - 6 



with the parameters rj s and £ s given by 

k 2z f _ a^puj 
ki z ki z 



(11) 



(12) 



where /io is the vacuum permeability. 

The transmission matrix for s and p polarizations 
across an interface has similar forms except for the sign 
of £ in the off-diagonal elements. Introducing a polariza- 
tion dependent parameter <; m , the transmission matrix 
for both polarizations can have an identical form 



Dx 



1 



Vm 



i 



Vn 



where m — (s,p) and q p = 1 and <; s 



Vm 

-1 



(13) 



B. Propagation matrix 

We now consider the propagation of light in a homo- 
geneous medium. It can be shown [33j that the electric 
or magnetic fields at z + Az can be related to those at 
the z position by a 2x2 propagation matrix 



P(Az) = 



-ik z Az 

o 







(14) 



C. Transfer matrix for multi- layer graphene 

For a stack of N graphene layers shown in figure [IJb) , 
the transfer matrix can be obtained by transmission ma- 
trices across different interfaces and propagation matrices 
in different homogeneous dielectric media. Denote the 
field coefficients on the left side of the leftmost graphene 
layer by a\ and b\ and those on the right side of the right- 
most graphene layer by a^+x and 6jv+i- The two sets of 
field coefficients are then related by a 2 x 2 transfer matrix 
Ai, namely 



ax 


= M 




> 




bN+i 



(15) 



with 



M = Dx^ 2 P(dx, 2 )D 2 ^ 3 P(d 2 , 3 ) ■ ■ ■ P(d N _x,N)D N ^ N+1 . 

(16) 
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D. Optical spectrum calculations 

With the transfer matrix, we can easily calculate the 
optical spectra such as reflection, transmission and ab- 
sorption for multi-layer graphene. Suppose that light is 
incident from left upon the multi-layer graphene with the 
reflection and transmission coefficients denoted respec- 
tively by r and t. It can be shown that these coefficients 
are given by the elements of M. 



t 



Mix' 
1 



Mil 



(17) 
(18) 



And reflectance and transmittance can be calculated for 
both s and p polarizations as 



Rs, P = \r 



s,p\ J 



(19) 
(20) 



where 



Vs = k(N+i)z/kiz, rjp = £ik(N+i)z/£N+ikiz- (21) 
Absorbance can then be readily obtained from 

A = 1 - R - T. (22) 

E. Optical conductivity of graphene 

For illustration and simplicity, in the present work we 
only consider the situation where the chemical potential 
of graphene \x is much larger than the temperature. In 
this situation, within the random-phase approximation 
the optical conductivity of graphene a(ui) is given by (l7r - 



EqC ll 



- 2) + - In 

7T 



n-2 



(23) 



Here, Q = tku/fi is the dimensionless frequency, a = 
e 2 /47reo^c (~1/137) is the fine structure constant and 
§{x) is the Heaviside step function. The first and sec- 
ond terms on the right side stem from the intraband and 
interband contributions, respectively. 



III. REFLECTION, TRANSMISSION AND 
ABSORPTION 

A. Single-layer graphene 

For a single graphene layer surrounded by two di- 
electrics of dielectric constants £\ and £2, suppose that 
light is incident from the dielectric medium of £\. The 
transfer matrix is no other than the transmission ma- 
trix across the interface, given by equations © and (|lip 




0.0 0.2 0.4 0.6 0.8 1.0 
Incident angle (ti/2) 



0.00 l 



0.0 0.2 0.4 0.6 0.8 1.0 
Incident angle (nil) 



FIG. 2. Reflectance, transmittance and absorbance of single- 
layer graphene as a function of the incident angle for Q > 2. 
(a) Reflectance and transmittance for s and p polarizations 
with £1 = 1 and £2 = 2.25. Dashed lines represent those 
in the absence of the graphene layer, (b) Absorbance for (i) 
ei = £2 = 1 and (ii) £1 = 1 and £2 = 2.25. (c) Same as (a) 
but for £1 = 2.25 and £2 = 1. The corresponding absorbance 
is given in (d). 



for p and s polarizations, respectively. From equations 
(HHJ) and (|2H)) . the reflectance and transmittance can be 
obtained as 



y/Ei COS 61 — ^/S2 cos 6*2 — a 



v /£l cos 6\ + y/E2~ cos9 2 +cr 

R yga/ cos 6> 2 - y/el/ cos 61 + a 2 
P ~ v/ei/cos02 + VeT/cosfli + CT ' 
T = 4 % /e 1 £ 2 cos 6<i cos6>2 , 2g . 
eTcos#i + costf 2 + 5f| 2 
4y£i£ 2 / (cosfli cos 6 2 ) , 

£2/ COS #2 + y/Slj COS 01 + o\ 



where 61 and 62 are incident and refracted angles, respec- 
tively and a = (t/bqc. By neglecting the higher-order 
terms of a, the absorbance is given by 



A. = 



A, 



cos 9iRe(cr) 



£l COS Vl + ^JS2 cos t/2 + o) 

4y / £iRe(5 : ) / cos Q\ 
feij cos 6 2 + \fz\l cos 9 1 + a) ' 



(28) 
(29) 



Obviously, for a single free-standing graphene layer under 
normal incidence its absorbance for fi > 2 is given by 
ira/ (1 + ira/2) 2 ~ ira. 

In figure [2j the reflectance, transmittance and ab- 
sorbance at different incident angles for single graphene 
are shown. Light is incident from the dielectric of E\. Re- 
flection and transmission are altered somewhat with re- 
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spect to the case without the graphene layer. For e\ < £2, 
the absorbance decreases monotonically with increasing 
incident angle for p polarization, while for s polarization 
it increases up to a maximum and then decreases mono- 
tonically. For £1 = £2, the absorbance for s polarization 
takes a universal value at normal incidence, about na 
and the maximal absorbance is 0.5 at an incident angle 
very close to tt/2. 

For £1 > £2, total internal reflection is expected. With 
the presence of the graphene layer, there is nearly no 
change in the critical angle. However, above the critical 
angle reflectance is no longer total (smaller than one). 
For s polarization, the absorbance around the critical 
angle is several times larger than the universal value of 
it a. This implies that the configuration of total internal 
reflection could be exploited in measurements of optical 
conductivity of graphene since it can suppress the signal- 
to-noise ratio considerably. 



B. Double-layer graphene 

Consider two graphene layers which are separated by 
a dielectric of £2 with a thickness of d. The dielectric 
constant of the leftmost dielectric is £1 and that of the 
rightmost dielectric is £3. The transfer matrix of the 
structure can be obtained from equation (If 61) as 



M = D 1 ^ 2 P(d)D 2 ^ 3 , 
whose elements are given by 
cos k 2z d 



M„ 



-A,. 



i sin k 2z d „ 



(30) 



(31) 
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FIG. 3. Reflectance, transmittance and absorbance of 
double-layer graphene as a function of the incident angle for 
Q, > 2. Two graphene layers have the same optical conduc- 
tivity, (a) Reflectance and transmittance for s and p polar- 
izations with d — O.lhc/fj, (corresponding to 0.132 /im for a 
typical chemical potential fi = 0.15 eV). Dashed lines repre- 
sent those in the absence of the graphene layers. The corre- 
sponding absorbance is given in (b). (c) Same as (a) but for 
d = 8hc/fi (corresponding to 10.5 fim for fi — 0.15 eV). The 
corresponding absorbance is given in (d). 



C. Multi-layer graphene 



where fi = 1, 2 and v — 1,2. The parameters A^m and 
are given by 



A^u.m — 1 + ( — l)** VrnVm 

B^, m = (-lTr hn + (r/ m + ^ m v'm) 

( 1)^ V?n£,m ^m£,m> 



where 



T)s = hz/hz, Vp = £lk2z/S2klz: 

£ s = anQw/ku, £,p = <Tk 2 z/saS20J, 

if s = hz/k 2z , r)' P = £2k 3z /£3k 2z . 

C s = o-^ oj/k 2z , ^ = ak 3z /e e 3 0J. 



(32) 
(33) 

(34) 



In figure [3j the reflectance, transmittance and ab- 
sorbance of double-layer graphene for e± = £3 = 1 and 
£2 = 2.25 are shown. For small d, despite different val- 
ues, the dependence of the reflectance, transmittance and 
absorbance on the incident angle is, in general, similar 
to that of single-layer graphene. At normal incidence, 
the absorbance is nearly twice the universal value of ira. 
For large d, however, there are oscillations in reflectance, 
transmittance and absorbance, which originate from the 
thin-film interference. 



For multi-layer graphene, the transfer matrix can be 
obtained from equation (|16[) . Reflectance, transmittance 
and absorbance can be then obtained from equations 
(HH), (J2QJ and ([22]) respectively. 

In figure [U the reflectance and absorbance of multi- 
layer graphene at normal incidence are shown. The struc- 
ture consists of identical graphene layers in air, separated 
equally by a distance d. At low frequencies the reflectance 
is close to one up to a certain cutoff frequency. Above the 
cutoff frequency, for small d the reflectance oscillates and 
drops rapidly to zero. For large d, however, there appear 
sharp reflection peaks for frequency above the cutoff fre- 
quency, resulting from the multiple interference by the 
graphene layers. 

The structure has zero absorbance for fl < 2 and shows 
remarkable absorption for il > 2. For small d, the ab- 
sorbance is nearly a constant, being 0.5 for ft > 2. With 
increasing d, multiple interference by graphene layers 
may play an important role, leading to sharp absorption 
dips. These interesting properties imply that multi-layer 
graphene has the potential as dark materials for achiev- 
ing lower reflection coatings and enhanced photodetec- 
tion [H. 




(a) 



(b) 



FIG. 4. Reflectance (solid line) and absorbance (dashed line) 
spectra of multi-layer graphene at normal incidence. There 
are totally 30 identical graphene layers, (a) d — O.lftc/ p 
and (b) d — 4hc/ fi (corresponding to 0.132 and 5.26 /jm for 
H — 0.15 eV, respectively). 



IV. PHOTONIC BAND STRUCTURE OF 
PERIODICAL GRAPHENE LAYERS 



When dielectrics are arranged in a periodical way 
to form so-called photonic crystals pfj| - l42| . electromag- 
netic waves should be strongly modulated by Bragg scat- 
terings, showing photonic band structures with well- 
defined photonic bands and even photonic bandgaps. For 
graphene layers stacking in a periodical way, photonic 
band structures should also be expected due to the in- 
troduced periodical modulations. 

The first structure considered is shown schematically 
in figure 03a). It is composed of identical graphene layers 
embedded into the interfaces of a one-dimensional pho- 
tonic crystal consisting of two dielectrics of e\ and e 2 
with a thickness of d\ and d 2 , respectively. For such a 
structure, its transfer matrix after propagating over one 
unit cell reads 



Mm = D 1 ^ m P(d 2 )D 2 ^ 1 , m P(d 1 ). 



(35) 



The photonic band structure can be then obtained 
from the diagonal elements of the transfer matrix [34j . 
cos(qd) — (Mn t7n + M22, m ) /2, where q is the Bloch 
wave-vector and d = d\ + d 2 is the lattice constant. It 
can be explicitly written as 



cos(qd) = cos(fci z di) cos(k 2z d 2 ) - -{r) m + rjj-) 
x sin(fci z di)sin(/c2 Z d 2 ) - A m , 



(36) 
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FIG. 5. (a) Schematic of a structure consisting of identical 
graphene layers placed at the interfaces of a one-dimensional 
photonic crystal, (b) Schematic of a structure consisting of 
identical graphene layers separated equally in air. (c) Pho- 
tonic band structure (solid lines) for the structure in (a). The 
parameters are ei = 1, £2 = 2.25, di = 3hc/n and d-z — 2hc//j, 
(corresponding to 2.63 and 3.95 fim for /1 = 0.1 eV, respec- 
tively). Dashed lines are the results in the absence of graphene 
layers, (d) Same as (c) but for the structure in (b). The pa- 
rameters are e — 1 and d = 4hc//j,. Dashed lines are simply 
the folded dispersion in air cj = qc. 



where 



A m = i£ m [(l + rim 1 ) sin(fci z di + fc 22 d 2 ) 
+Cm (l - Vm 1 ) sin(k lz di - k 2 zd 2 )} 

tm 

+5— sin(fci z di) sin(fe 2 d 2 ). 



(37) 



Without the term A m , equation (|36| reduces to the pho- 
tonic band structure of a one-dimensional photonic crys- 
tal [13. 

For identical graphene layers separated equally by a 
distance d in air as shown in figure EJb), the transfer 
matrix simply reads Ai m = D m P{d), where 



D„ 



1 ~\~ £m/2 <vm£m/2 



(38) 



The photonic band structure of the structure is then 
given by 

cos(qd) — cos(k z d) — sin(fc z d), (39) 

which is identical to that given in reference [43[ . 

In figure 03 the photonic band structures of period- 
ical graphene layers are shown for the propagation di- 
rection perpendicular to the graphene layers. For both 
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structures the first photonic band starts from a certain 
nonzero cutoff frequency, different from conventional di- 
electric photonic crystals. This cutoff frequency corre- 
sponds exactly to that observed in the reflection spectra 
shown in figure [4] Below the cutoff frequency, high re- 
flection or low transmission is expected, which was also 
observed in numerical simulations of a stack of graphene 
layers separated by dielectric slabs [44[ ■ It is known that 
in the low-frequency limit periodical metallic structures 
can be considered as bulk metals with a depressed effec- 
tive plasma frequency [45|,[46[. Thus, periodical graphene 
layers can also be regarded as a bulk metal with an ex- 
tremely low effective plasma frequency. 

For the structure shown in figure EJa), the photonic 
crystal in the absence of graphene layers displays well- 
defined photonic bands and bandgaps. With the in- 
troduction of graphene layers, however, both photonic 
bands and bandgaps are modified. For example, for the 
photonic crystal in the absence of graphene layers, there 
should be no bandgap between the second and third pho- 
tonic bands since this photonic crystal is a quarter-wave 
stack. In the presence of graphene layers, however, a mini 
photonic bandgap opens up. 

For the second structure shown in figure EJb) , there 
should be no photonic bandgaps in the absence of 
graphene layers. In the presence of graphene layers, how- 
ever, a series of mini photonic bandgaps appear owing to 
the multiple interference by graphene layers. It is known 
that for frequency within photonic bandgaps light prop- 
agation is forbidden poMi^ . For a structure consisting 
of finite graphene layers, this will cause strong reflection 
for frequency located into the mini photonic bandgaps, 
as can be clearly seen from figure EJb). 



V. PLASMONS 

A flat metal surface can support surface plasmons 
[24Tj26j which are transverse magnetic (TM) electromag- 
netic waves coupled with collective oscillations of surface 
charges. Surface plasmons can propagate along the metal 
surface with the fields decaying exponentially away from 
the both sides of the surface. In doped or gated graphene, 
free carriers can also support plasmons [17H22| . Owing to 
the two-dimensional nature and unique electronic band 
structure, graphene can support not only TM but also 
transverse electric (TE) plasmons The later does 

no exist in conventional metal surfaces. 

For a graphene layer surrounded by two dielectrics as 
shown in figure [IJa) , the transfer matrix is simply the 
transmission matrix. From equation (|17[) . the reflection 
coefficient of the system can be obtained. The condi- 
tion for the existence of plasmons is that the reflection 
coefficient has poles, namely 



1 + Vm + ?ro£m = 0, 



(40) 



where the subscript m stands for s and p polarization, 
corresponding to TE and TM modes, respectively. 



From equation (|40)) . the dispersion of TM plasmons 
can be obtained as 



si 



S2 



^Q 2 - £l W VO 2 - ^ 2 



r(Q)/goc 

n 



(41) 



with Q = hcq/fi, where q (= k x ) is the wave- vector of 
the plasmons. Obviously, TM plasmons can exist if the 
imaginary part of a is positive. From the above equation, 
TM plasmons are far below the light line, i.e., q 3> co/c. 
Thus, in this non-retarded regime, the dispersion of TM 
plasmons is simplified to 



Q 



(si + e 2 ) Cl(s c)/ia. 



(42) 



For SI > 2, a has a real value as well, leading to a strong 
loss due to interband excitations. For small q, the dis- 
persion of TM plasmons reduces to 



n 



(ei + £2) 



Q. 



(43) 



which shows the known ^-dependence [13, 13 • 

From equation (|40| , the dispersion of TE plasmons is 
given by 



VQ 2 - £i^ 2 + VQ 2 - = —SI, 

e c 



(44) 



which is the same as that given in reference |l9j. TE 
plasmons can exist if ei = £2 and the imaginary part 
of a is negative (for SI > 1.667). Note that the term 
on the right side of the above equation is very small in 
the frequency window 1.667 < SI < 2. As a result, the 
dispersion of TE plasmons should be below but very close 
to the light line u — qcj *Je\. 

We now consider a structure where a graphene layer 
is apart from a dielectric substrate of £2 over a distance 
of d, as shown schematically in the inset of figure [6] For 
the structure, TE plasmons do not exist for finite d and 
thus we only discuss TM plasmons. For p polarization, 
the transfer matrix of the structure can be obtained from 
equation (|16|) . namely 



where 



D 



D' 



P 



M = DP(d)D', 



1 - v P + 1 + vp - iv 

1 + nLl 



-iki z d q 

e lkl * d 



"Hp 



(45) 

(46) 
(47) 
(48) 



with 



£ P = (rk lz /e siuj, rj' p = e 1 k 2z /e 2 ki z . (49) 
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1000 



FIG. 6. Dispersion of TM plasmons for a structure shown 
in the inset where the graphene layer is apart from a dielec- 
tric substrate of £2 = 10 over a distance of d and the other 
dielectrics are air with ej = 1. Blue, green and red lines are 
the results for d — 0.005, 0.01 and 0.1, respectively. The dot- 
ted line represents the dispersion for a free-standing graphene 
layer in air and the dashed line corresponds to the case d — 0. 



With the transfer matrix, it is easy to obtain the reflec- 
tion coefficient from equation (jTTf . The condition for the 
existence of plasmons thus reads 



(1 



-ik\ z d 



(1 - rQ e^ d = 0. (50) 



Since the dispersion of graphene plasmons lies far below 
the light line, the non-retarded condition q 3> w/c still 
holds, leading to rj' p ~ £i/£2 and £ P — icrq/eoSiUJ. Thus, 
plasmon dispersion can be simplified as 



2ei (ei 



■ £2 



(ei + £2) + (£1 - £2) e 



-2Qci 



EqC fl ' 



(51) 



where d — d^i/Tic. 

In figure [SJ the dispersion of TM plasmons for the 
structure shown in the inset is given. The separation 
of the graphene layer from the dielectric substrate of £2 
affects considerably the dispersion. For small q, the dis- 
persion takes that for the case of d = 0. For large q, it 
approaches that for the free-standing case. This disper- 
sion interchange with increasing q can be understood by 
the fact that the fields of plasmons decay exponentially 
into the surrounding media, as schematically depicted in 
the inset. For small q, the decay length is much larger 
than d such that the fields concentrate dominantly in the 
substrate. As a result, the dispersion should take that 
for the case d ~ 0. For large q, the fields decay very 
fast such that the decay length is much smaller than d. 
In this situation, the fields could not sense the substrate 
dielectric layer. 

Obviously, the dispersion interchange can be tuned by 
d. The change from one kind of dispersion to the other 



occurs faster for large d than for small d. Practically, we 
may adopt this dispersion interchange to tune the disper- 
sion of TM plasmons, or in other words, the refractive 
index of the plasmons. From equation (jl!?)) , the corre- 
sponding refractive index of TM plasmons for a graphene 
layer surrounded by two dielectrics of £\ and £2 is given 
by 



qc/u 



£1 +£2 
ia/£oC 



(52) 



With the structure shown in the inset of figure |6l the 
refractive index of TM plasmons can thus be tuned from 
— 2ei/ (ia-/£o c ) t° — (£1 + £2) / («cr/£o c )- This offers prac- 
tically a simple approach to manipulate the dispersion of 
TM plasmon or the refractive index by changing d. 



VI. CONCLUSIONS 

In the present paper, we developed a transfer ma- 
trix method for optical calculations in non-interacting 
graphene layers. With the framework of this method, 
the transfer matrix for various graphene layers can be 
obtained. With the transfer matrix, reflectance, trans- 
mittance and absorbance spectra of graphene layers can 
be easily obtained. In addition, photonic band struc- 
tures for periodical graphene layers and even graphene 
plasmons can be studied in a rather simple way. 

Using the transfer matrix method, we studied the opti- 
cal properties such as reflection, transmission and absorp- 
tion for single-, double- and multi-layer graphene. We 
showed that the configuration of total internal reflection 
in single-layer graphene and thin-film interference effects 
in double-layer graphene could be exploited to enhance 
the light absorption. For multi-layer graphene, there ex- 
ists a cutoff frequency below which the reflectance is as 
high as one. For a small spacing distance, the absorbance 
is as high as 50% for ft > 2. With increasing spacing dis- 
tance, sharp reflection peaks and absorption dips appear 
owing to the multiple interference by graphene layers. 

We applied the transfer matrix method to structures 
consisting of periodically arranged identical graphene 
layers. The structures are characterized by photonic 
band structures with well-defined photonic bands and 
bandgaps. We revealed that these structures can be re- 
garded as bulk metals with an extremely low effective 
plasma frequency. 

Finally, we discussed plasmons in a graphene layer 
which is apart from a dielectric substrate. We found that 
the plasmon dispersion can be tuned by the separating 
distance between the graphene layer and the dielectric 
substrate. Our results show that the transfer matrix 
method could serve as a versatile tool to study optical 
properties in graphene layers. 
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